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Abstract: A strategy for computing the ^ip anomalous dimension at the fixed point in infrared- 
conformal gauge theories from lattice simulations is discussed. The method is based on the scaling 
of the spectral density of the Dirac operator or rather its integral, the mode number. It is relatively 
cheap, mainly for two reasons: (a) the mode number can be determined with quite high accuracy, 
and (b) the il^ip anomalous dimension is extracted from a fit of several observables on the same 
set of configurations (no scaling in the Lagrangian parameters is needed) . As an example the ip^ 
anomalous dimension has been computed in the SU(2) theory with 2 Dirac fermions in the adjoint 
representation of the gauge group, and has been found to be 7* = 0.371(20). In this particular case, 
the proposed strategy has proved to be very robust and effective. 
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1 Introduction 

Non-Abelian gauge theories with minimally coupled Dirac fermions generally produce spontaneous 
chiral symmetry breaking if the number of fermions (or flavors) is not too large. When the number 
of flavors is increased,^ chiral symmetry is restored and the beta function develops a non-Gaussian 
infrared fixed point [1], indicating that an approximate dilatation symmetry is realized at large 
distances. Gauge theories in this phase are usually referred to as infrared- conformal gauge theories 
(IR-CGTs). Deformed (with relevant operators) IR-CGTs have been proposed in the past years as 
interesting models for strongly-coupled physics beyond the Standard Model, for instance in walking 
or conformal technicolor scenarios (see [2-11], and references therein). The anomalous dimension 
7* of the fermion mass operator ip^p at the IR-fixed point plays a very important role in all those 
models. Being constrained to be in between (as in the perturbative Bank-Zaks fixed point [1]) 
and 2 (from the unitarity bound [12, 13]), it is required to be about 1 for interesting models. In 
this regime the IR-CGT is strongly coupled and a determination of 7^ from first principles can be 
obtained only by means of lattice simulations. 

The theory considered in this paper is the SU(2) gauge theory with 2 Dirac fermions in the 
adjoint representation of the gauge group (SU(2) -I- 2adj). It is IR-conformal as confirmed by a 
relatively large literature and several complementary analysis strategies [14-34] . Although it is not 
clear whether this theory will be useful for building realistic technicolor models, it represents the 
ideal playground for testing new analysis methods. The ipip anomalous dimension has been already 
estimated using different techniques. The first very rough estimate 0.05 < 7* < 0.20 was published 
in [20] using finite-size scaling of mesonic observables (similar to the one given in [35]). A somewhat 
larger value 7* = 0.22(6) was found in [19], fitting a power of the fermion mass to the string tension. 
These first estimates must be taken with a grain of salt, since the systematic errors were still not well 
understood at that time. Using the Schrodinger Functional renormalization scheme, the authors 

^If the number of flavours is large enough, asymptotic freedom is destroyed. Only asymptotically free theories 
will be considered here. 
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of [30] quote 0.31(6) and the authors of [29] quote 0.05 < 7* < 0.56. As discussed in [29] large 
systematic errors are expected with this technique, because of the difficulty in localizing the IR-fixed 
point. Monte Carlo Rcnormalization Methods [32] have similar problems leading to a very wide 
allowed range —0.6 < 7* < 0.6. A more sophisticated analysis of the finite-size scaling of mesonic 
masses and decay constants [34] yields 7* = 0.51(16). The central value is larger than any other 
determination, but the error is also quite large (about 30%). 

The method proposed in this paper to measure 7* is based on the observation that the spectral 
density p of the Dirac operator is a power of the eigenvalue w . Its exponent is related to the 
anomalous dimension through the relation [36-38] 

p(w) poAi^wTT^ + . . . , (1.1) 

at the leading order for small eigenvalues (/i is the renormalization scale and is a dimensionless 
constant). When scale invariance is broken by a small mass for fermions (mass-deformed IR- 
conformal gauge theory, IR-mCGT), eq. (1.1) is valid only in an intermediate range of eigenvalues, 
that is not known a priori and must be determined empirically. 

In principle this method is limited by the need to have a reliable infinite-volume extrapolation, 
and to be close enough to the chiral limit. However we will see that the fermion mass does not need 
to be too small in order to see a power law in the spectral density for a wide range of eigenvalues. 
This surprising empirical observation, combined with the fact that the spectral density (or rather 
its integral) can be computed with high accuracy, will lead to a quite precise determination of the 
anomalous dimension, that is anticipated to be 7* = 0.371(20). Although this result is obtained for 
a particular theory, the method can be exported with no modifications to any other gauge group 
or matter content. A precursor of the presented strategy can be found in [38]. While this work 
was being finalized, an analysis of the Dirac spectral density in the SU(3) theory with 12 fermions 
in the fundamental representation appeared in [39] , although very few eigenvalues were considered 
there. In this work the 50 lowest eigenvalues have been discarded on a 64 x 24^ lattice, and about 
2000 eigenvalues have been considered in order to extract the ij^tp anomalous dimension. 



2 Spectral density 

The validity of eq. (1.1) is discussed in this introductory section. The eigenvalue density of the 
massless Dirac operator If) (or Dirac spectral density) is defined as 

p(c.)= hm i^(<5(c.-c.,.[A])) , (2.1) 



k 



where ja;fc[A] is the fc-th eigenvalue of Ip at fixed gauge configuration, in a finite box of volume 
V. At the root of eq. (1.1) is the fact proved in [40], that both the eigenvalue w and the Dirac 
spectral density p{uj) renormalize multiplicatively. If 7(17) is the ipil) anomalous dimension at a given 
coupling (7, the Dirac spectral density p{ijj) has anomalous dimension 7(5), while the eigenvalue u 
has anomalous dimension —7(g). 

An IR-CGT is characterised by the existence of an IR-fixed point for the running coupling 
g = g*- The leading behaviour in eq. (1.1) is obtained from the renormalization-group equation 
exactly at the fixed point [36-38]: 



d 

(1 + 7*)^^ - 3 -f 7* 



piuj) = , (2.2) 



which depends only on the ^pip anomalous dimension 7^ = 7(5*) at the IR-fixed point. However IR- 
CGTs are asymptotically free at high energies and reach the IR-fixed point only at asymptotically 
low energies. Hence eq. (2.2) and its solution (1.1) are valid only for low eigenvalues. 
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A mass term for fermions destroys scale invariance. If the mass is small enough, the renormal- 
ization flow starts from the vicinity of the Gaussian fixed point at very high energy, moves towards 
the IR-fixed point, stays in its vicinity for some time, and is eventually driven away from it by the 
mass term. A smaller mass corresponds to a longer time spent in the vicinity of the IR-fixed point. 
As long as the renormalization flow stays close to the IR-flxed point, approximate scale invariance 
is generated. 

Therefore in IR-mCGTs one can identify three regions for the eigenvalue density: 

1. An intermediate region of eigenvalues Ojn < uj < Cjuv- This region is dominated by 
the vicinity of the IR-flxed point. The power law (l.f) is valid only in this intermediate range 
of eigenvalues. 

2. A region of high eigenvalues lj > Coijv- Around uj ~ ojuv the spectral density begins 
its transition to the region dominated by asymptotic freedom, which will eventually lead to 
its asymptotic form p{uj) oa uj^. At these large eigenvalues, the spectral density is essentially 
insensitive to the fcrmion mass, hence Ljjjv has a well-deflned nonvanishing chiral limit. 

3. A region of low eigenvalues < a; < u)ir. The spectral density feels the fermion mass 
and is determined by the details of the dynamics of the IR-mCGT. No analytical model is 
available here a priori. This region disappears in the chiral limit {uiir goes to zero). Since 
at these low energies the physics decouples from the UV scale determining the transition to 
the asymptotically free regime, and since eigenvalues have the same scaling dimension of the 
mass, the IR scale ujm must be proportional to m: 

UJIR (X m (2.3) 

at the leading order in the chiral limit. 



3 Method 

The proposed method to extract the ipip anomalous dimension from lattice simulations is described 
in this section. Although the Euclidean Dirac operator is diagonalizable in the continuum theory, it 
is not generally so on the lattice, and one can more conveniently study the positive-definite operator: 

M = {p + m)^p + m) = - . (3.1) 

Following [40], the mode number per unit volume P(ri) is defined as the number of eigenvalues of 
M lower than fi^ divided by the volume: 

= 2 p{oj) dw . (3.2) 
Jo 

In an intermediate region oj'jj^ + < il? < oJijv + the mode number per unit volume is: 

D{n) = 2 / p{uj) duj + 2 p{uj) duj . (3.3) 

The first term will be named ^^{m). It is an unknown additive constant, independent of J7, that 
depends only on the mass m (remember that uim oc m). As discussed in sec. 2, in the second 
integral the spectral density can be approximatively replaced by the power law in eq. (1.1). Putting 
it all together: 

D{n)~Do{m) + ^{l + j,)poP^{n^-m^)^ . (3.4) 
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Notice that the prefactor and the exponent of the power law do not depend on the fermion mass, 
while the additive constant Pq does (and, in particular, vanishes in the chiral limit). Since the mode 
number per unit volume is renormalization-group invariant [40] and has the dimension of a mass to 
the 4th power, the additive constant must scale like 

Doim) cx Mis , (3.5) 

where Mps is the mass of the isotriplet pseudoscalar meson (or any other particle mass). 

All formulae written so far are expressed in terms of renormalized quantities. In order to 
compare eq. (3.4) to lattice data, one has to trade the renormalized quantities with the bare ones. 
This can be done by replacing (see [40] for details) 

" ^ ^ , m-^ — ^ —mpcAC , (3.6) 

Zj p Zjp Zi p 

where Za and Zp are the renormalization constants of the isovector axial current and the 
isovector pseudoscalar density P° respectively, and mpcAC is the quark mass as defined from the 
hare PCAC (partially-conserved axial current) relation 

d^^Al = mpcAcP" . (3.7) 

The mode number does not need to be redefined, since it is renormalization-group invariant. Putting 
everything together, the mode number per unit volume in lattice units as a function of the hare 
eigenvalue O becomes 

a-^v{^) ~ a-^VQ + A\(aVLf ~ {amf\ ^ , (3.8) 

for a suitable definition of the dimensionless constant A. The parameter m appearing in the previous 
formula is related to the hare PCAC mass through the relation 

m = ZA-mpcAC ■ (3-9) 

The proposed strategy consists in computing the mode number per unit volume hy means of 
lattice simulations, and in extracting the iptp anomalous dimension hy fitting eq. (3.8) to the lattice 
data. Since the validity range ^ir < f2 < fluv of eq. (3.8) is not known a priori, it has to he 
determined hy studying quality and stahility of the fit procedure. 

The mode number per unit volume can be computed essentially in two ways. 

1. The eigenvalues of AI can be explicitly computed starting from the lowest one using the 
Chebyshev accelerated subspace iteration method described in detail in [41] . Since the number 
of eigenvalues below some fixed value fi^ grows linearly with the volume, this method is 
unpractical for large volumes. 

2. Alternatively (following ref. [40]) one can define the projector P{il) over the eigenspaces of 
M corresponding to lower eigenvalues than fl^ , in terms of which the mode number per unit 
volume is 

D{n)^ lim -(trP(r2)). (3.10) 
The projector can be approximated by a suitable rational function: 

p{n)c,h{x)\ x^i-j^^, (3.11) 

where h{x) is a polynomial, and fl^ is a parameter of order 57 defined in eq. (A. 5). The 
trace in eq. (3.10) can be estimated stochastically. The error due to the approximation in 
eq. (3.11) can be estimated a posteriori, once the spectral density is reconstructed by means 
of the fit procedure (see appendix A). This technique is more effective for larger volumes, and 
is described in detail in [40]. 
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Set 


lattice 


airiQ 


ampcAC 


aMpg 


References 


SI 


64 X 


-1.15 


0.11844(63) 


0.6414(41) 


[23] 


S2 


64 X 32^ 


-1.15 


0.11790(36) 


0.6386(15) 


[23] 


S3 


64 X 243 


-1.18 


0.05507(89) 


0.3374(63) 


[21] 



Table 1. Summary of the three sets, used for computing the mode number. The PCAC and isotriplet 
pseudoscalar meson masses have been published in the cited papers, and are given here just for the reader's 
convenience. 

4 Analysis and results 

In order to illustrate the proposed technique, let us consider the SU(2) + 2adj theory, which is in 
the conformal window. 

Numerical simulations have been performed at fixed = 2.25 with Wilson fermions. I have 
used three sets (table 1) of configurations generated with the HiRep code [16] and already used in 
previous works to compute the mass spectrum (mesons, glueballs and string tension) [18-21, 23]. 
The first set (SI) corresponds to bare mass arriQ ~ —1.15 on a 64 x 24"^ lattice, and includes 25 well- 
decorrelated configurations. The second set (S2) corresponds to the same bare mass amo — —1.15 
but on the larger 64 x 32^ lattice, and includes 20 well-decorrelated configurations. The third 
set (S3) corresponds to a lighter bare mass arriQ = —1.18 on a 64 x 24^^ lattice and includes 20 
well-decorrelatcd configurations. 

The ^1}) anomalous dimension will be extracted by fitting eq. (3.8) to the mode number per unit 
volume obtained from sets SI and S3. Set S2 will be used to check that the finite-volume effects 
are under control at least for the higher mass. Moreover the stability of the fitting procedure will 
be checked against the introduction of a subleading term in eq. (3.8). 

Before moving to the details of the analysis, the results of the following subsections are antic- 
ipated and summarized in table 2 for the reader's convenience. The ip^ anomalous dimensions in 
table 2 are all compatible; however, the first one 7* = 0.371(20) obtained from set SI is quoted 
as the final results. This choice is dictated by the fact that the finite-volume effects have been ex- 
plicitly checked for set SI. Even considering the weighted average between the two determinations 
obtained in the same range 0.091 < afi < 0.18, one would get the very similar result 7* = 0.365(19). 



Subsection 


Data set 


Fitting function 


Fit range 


7* 


4.1 


SI 


(3.8) 


0.091 <an< 0.18 


0.371(20) 


4.2 


SI 


(4.2) 


0.091 <an< 0.29 


0.355(23) 


4.4 


S3 


(3.8) 


0.091 <an< 0.16 


0.364(92) 


4.4 


S3 


(3.8) 


0.091 <an< 0.18 


0.325(50) 



Table 2. Determinations of the anomalous dimension from: (a) set SI (amo = —1.15); (b) set SI 
(amo = —1.15) including a subleading correction; and (c),(d) set S3 (amo = —1.18) with two different 
fit ranges. The three determinations are compatible, in the sense that the la regions overlap. The first 
determination is quoted as the final result. 

4.1 Set SI: determination of 7* 

The lowest 200 eigenvalues of M are explicitly computed. This allows us to reconstruct the mode 
number per unit volume up to afJ = 0.1604. Above this eigenvalue, the projector method is used. 
Results are listed in table 3. 

For aVl ~ 10 the mode number per unit volume saturates at the value 120""* (12 is the dimension 
of the vector space spanned by the pseudofermions in a single site) . The saturation is a pure lattice 
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artifact, so only values of afl for which the mode number per unit volume is less than 0.12a~'* are 
considered (one hundredth of the total number of eigenvalues) . 

Eq. (3.8) is fitted to the lattice data. Although the PCAC mass is known, but no determination 
of Za is available, the parameter am which determines the mode-number gap can not be recon- 
structed from eq. (3.9) and is considered as a fitting parameter along with a~'*Poj A and 7*. From 
the discussion in sec. 3 it is also clear that the fitting function in eq. (3.8) can be used only in an 
intermediate range of eigenvalues, that must be determined by studying the stability and quality 
of the fit. 

Determination of the fit-range lower end. The fit works well with a x^/dof of order 1 in 
the quite large range 0.08 < aft < 0.4. Choosing the fit-range higher end to be 0.18, the lower end 
is systematically increased from 0.08 to 0.1 (fits S1:F1 to S1:F7 in table 6). The fit parameters are 
shown as functions of the fit-range lower end in fig. 1. Increasing the fit-range lower end, the fit 
parameters keep shifting up to afi// ~ 0.91 where they reach a plateau. The latter value is chosen 
as the fit-range lower end. Notice that at the lower end of the fit range P(0.092/a) > SDq, and the 
additive constant becomes soon negligible for larger eigenvalues. 




0.088 



E 0.084 

CO 




0.08 0.085 0.09 0.095 0.1 



-2e-05 



0.08 0.085 0.09 0.095 0.1 



Figure 1. Parameters in eq. (3.8), fitted to set SI in the range u^Il < fflfi < 0.18 for various values of the 
lower end aQ.L- While the fit range is reduced the fit parameters shift, and reach plateaux for afii — 0.091. 

Determination of the fit-range higher end. The fit-range higher end is also systematically 
investigated. Eq. (3.8) is fitted to the data in the range 0.091 < afl < aVln, where the higher end 
a^H is gradually lowered from 0.6 to 0.1 (fits S1:F8 to S1:F26 in table 6). Ahhough the x^/dof 
becomes less than 1 at about aVln — 0.4, the fit parameters keep shifting up to aVLn — 0.18 where 
they reach a plateau (see plots in fig. 2). The latter value is chosen as the fit-range higher end. 

Summarizing, the chosen fit range is 0.091 < aVl < 0.18 (fit S1:F4 in table 6), which yields a 
mode number per unit volume of the form (fig. 3) 

a-^i?{n) ~ 1.31(78) X 10"^ + 0.532(38) [(ar^)^ 0.0826(16)^] i+o-^tkso) . (4.1) 

The -ipip anomalous dimension is determined to be 7* = 0.371(20). Notice that the chosen range fit 
corresponds to about 2000 eigenvalues on the 64 x 24^^ lattice. 
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Figure 2. Parameters in eq. (3.8), fitted to set SI in the range 0.091 < aO, < aQn for various values of the 
liigher end ailn- While the fit range is reduced the fit parameters shift, and reach plateaux for aQn — 0.18. 



4.2 Set SI: corrections to the leading power 

It is interesting to check the stability of the fit result (4.1) against the introduction of a subleading 
power-law contribution in eq. (3.8), which becomes 



(4.2) 



Since eq. (4.1) describes very well the data in the range 0.091 < ail < 0.18, one needs to enlarge 
the fit range in order to discriminate the subleading contribution. Fitting the data to eq. (4.2) in 
the range 0.091 < afl < 0.29, one obtains 



a-^D{n) ~ -0.06(68) X lO""^ + 0.529(49) [(all)^ - 0.0798(14)^] 1+0.355(23) + 

+ .16(ll)[(al})2 - 0.0798(14)2]-926(90) 



(4.3) 



Notice that by including the subleading contribution, the estimate for the ipip anomalous dimension 
in this range has increased from 7* = 0.3212(42) to 7* = 0.355(23), becoming compatible with the 
determination of the previous subsection (in the sense that the la regions overlap). 
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2 2 1/2 

Figure 3. Mode number per unit volume for the set SI {arriQ = —1.15 on a 64 x 24"^ lattice): lattice 
data and fit result in log-log scale. The reference fit is S1:F4 in table 6. The parameters in the axis labels 
have been chosen to be a~'^VQ = 1.31 x 10~^ and am — 0.0826 (best-fit results). The black points are the 
data computed by numerical simulations. The red line is the best fit to eq. (3.8), while the orange band 
corresponds to the Icr region. The blue dashed lines delimit the data used for the fit. 

4.3 Set S2: finite-volume effects 

As analyzed in [23], meson masses computed on the set SI {arriQ = —1.15 on 64 x 24^) are identical 
to the ones computed on the set S2 (amo = —1.15 on 64 x 32'^), within the statistical errors that 
are of the order of 0.5%. It is reasonable to expect that finite-volume effects are under control 
also for the mode number. However this is explicitly checked by computing the mode number per 
unit volume using the projector method for few values of afl. The agreement is always within 
la as shown in table 4. Since larger finite-volume effects are expected for lower eigenvalues, we 
can conclude that the finite-volume effects for the set SI are always negligible with respect to the 
statistical errors for oJl > 0.086. 

4.4 Set S3: lighter mass 

The set S3 (amo = —1.18 on 64 x 24'^) is used to check the stability of the 'ipip anomalous dimension 
while going closer to the chiral limit. For this set no detailed investigation of finite- volume effects 
is available. However the isotriplet pseudoscalar meson is expected to be about 10% lighter than 
in infinite volume (see analysis in [23]). Similarly one has to expect sizable finite- volume effects 
also for the spectral density at low eigenvalues, while for larger eigenvalues the finite-volume effects 
become smaller. I will work under the assumption that the finite-volume effects are comparable in 
the two sets SI and S3 at fixed eigenvalue. Therefore the analysis is restricted to the safe range 
an > 0.086. 
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The mode number per unit volume for this set has been computed using the projector method, 
and results are listed in table 5. Eq. (3.8) has been fitted to the data in the range 0.091 < ail < ailn, 
for several values of the right end aflff (fits S3:F1 to S3:F5 in table 6). All fits have x^/dof of order 
1, and give compatible results (the la regions overlap). The value of the ijj^ anomalous dimension 
obtained from the largest fit range is 7* = 0.325(50) (fit S3:F5 in table 6). Both the amplitude A 
and the anomalous dimension 7*, that are mass- independent, are compatible with the ones obtained 
from the heavier mass (fit S1:F4 in table 6), in the sense that the la regions overlap. 

5 Conclusions 

The ■0V' anomalous dimension of the SU(2) gauge theory with 2 Dirac fermions in the adjoint 
representation (SU(2) + 2adj) is extracted from the mode number per unit volume of the operator 
M = m? — : which is related to the Dirac spectral density piuj) by 



The mode number per unit volume is computed by means of lattice simulations, using the methods 
described in [41] and [40]. 

In an IR-conformal gauge theory (IR-CGT), the spectral density at small eigenvalues follows a 
power law, and so does the mode number. The exponent is related to the V'^ anomalous dimension 
7* at the IR-fixed point. If conformality is explicitly broken by a mass for the fermions, the spectral 
density is expected to follow the power law only in an intermediate range of eigenvalues. In this 
intermediate region, the mode number per unit volume (in lattice units) has approximatively the 
form 



This work proves that it is possible to use the previous formula to extract the V'V' anomalous 
dimension from lattice simulations. The strength of this method relies on the following facts: 

• The mode number can be obtained with quite high accuracy through lattice simulations, even 
with few configurations (25 configurations have been used for the higher mass), since it is 
an extensive quantity. Its computation does not involve fitting procedures (like for particle 
masses) . 

• Since the finite-volume effects on the mode number are smaller for larger eigenvalues, useful 
information can be extracted also from not-so-large volumes. 

• The "0^ anomalous dimension is extracted from different observables (i.e. the mode number 
at different eigenvalues) computed on the same set of configurations. All the other known 
methods are based on the analysis of the scaling of a fixed observable with some parameters 
in the action (typically the mass, the volume or both). 

It is also surprising to observe that the region controlled by the power law is already quite wide 
at some intermediate mass {aMpg ~ 0.5). This last observation has an empirical nature and is in 
principle model-dependent. 

The -i/jV anomalous dimension of the SU(2) -f 2adj is found to be 7* = 0.371(20) at fixed 
P = 2.25. The stability of this result has been checked by lowering the fermion mass and by 
including a subleading power in the fitting function. An analysis of the continuum limit will be 
attempted in the future. It is worth recalling again that, even though the presented result is 
obtained for a particular theory, the method can be exported with no modifications to any other 
gauge group or matter content. 




(5.1) 



(5.2) 
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A Projector-approximation error estimate 

It is worth recalling a few essential facts about the procedure for the approximation of the projector 
P(r2) (for details the reader should refer to [40]). Let P{x) be the min-max polynomial of degree n 
that minimizes the deviation 



max |l-y^P(y)| , (A.l) 

e<y<l 



and define the function 



h{x)^l[l-xPix')] . (A.2) 

The polynomial P{x) approximates the function x^^^^ in the range e < x < 1, whilst h{x) approxi- 
mates the function 0{—x) in the range -^/e < \x\ < 1. Two degrees of approximation have been used 
in this paper (depending on the required precision): (a) a polynomial of degree 32, with e = 10~^ 
and 6 ~ 4.35 x lO"'^; (b) a polynomial of degree 100, with e = 10"^ and 6 ~ 5.20 x 10"''. 

The error due to the approximation of the projector P(51) with the rational function /i(X)'* has 
the following spectral representation: 

duj [e{fl-u})-h{x^)]a-^p'iu}) , (A.3) 
= (A.4) 

The quantity $7* is defined as 



dx ^T^Hx^ ■ (A.5) 



The error A in eq. (A.3) has been estimated by using the functional form of the mode number 
per unit volume obtained from the best fit: 

_4_, , fa-4pQ + o.532[(al7)2 - 0.08262]Ti7T for amo = -1-15 , 
a 1^(11) ~ < , . (A. 6) 

[a-^Da + 0m7[{an)^ - 0.0482]rai5 for amo = -1.18 

In all cases the error A has been checked to be smaller that the statistical error on a^'^D{il). 
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5.282(13) X 10-3 


0.2411 


7.217(16) X 10-3 


0.2641 


9.810(19) X 10-3 
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9.5130(62) X 10-2 
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Table 3. Set SI (64 x 24^ /3 = 2.25 amo = —1.15). Mode number per unit volume, computed from the 
eigenvalues for af2 < 0.11 and with the projector method for afi > 0.11. 



an 


a-4p(f]) @ 64 X 243 (SI) 


a-4p(Sl) @ 64 X 323 (^2) 


0.086 


2.05(12) X 10-5 


1.974(65) X 10-5 


0.092 


5.91(21) X 10-5 


5.90(23) X 10-5 


0.098 


1.1122(28) X 10-4 


1.1105(26) X 10-4 


0.104 


1.811(31) X 10-4 


1.827(19) X 10-4 



Table 4. The mode number for the set S2 (64 x 32^ /? = 2.25 amo = —1.15), computed with the projector 

method, is compared to the mode number for selected points of the set SI (64 x 24^ /3 = 2.25 amo = -1-15). 
For all the considered eigenvalues, the finite-volume effects are negligible. 
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1.5146(77) X 10-3 
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1.896(11) X 10-3 
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2.337(12) X 10-3 


0.180 


3.423(16) X 10-3 



Table 5. Set S3 (64 x 24^ ^ = 2.25 amo = —1.18). Mode number per unit volume, computed with the 
projector method. 
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0.449(68) 


8.50(21) X 10"^ 


0.415(37) 


17 


0.22 


S1:F16 


0.091 < afi < 0.16 


0.13(10) X 10"" 


0.533(59) 


8.25(22) X 10"^ 


0.370(29) 


19 


0.33 


S1:F17 


0.091 < an< 0.17 


1.26(83) X lO"'^ 
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20 
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22 


0.33 


S1:F19 


0.091 <an< 0.21 


0.38(72) X 10"^ 
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25 
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0.641(13) 
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26 
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27 
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30 
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0.091 < an < 0.5 
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33 


2.4 
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0.091 < an < 0.6 


3.58(52) X 10"^ 


0.5805(13) 


8.54(83) X 10"^ 
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35 


10 


S3:F1 


0.091 < an < 0.13 


-0.04(16) X 10"^ 


0.85(31) 


0.12(45) X 10"^ 


0.25(18) 


6 


1.4 


S3:F2 


0.091 < an < 0.14 


0.02(11) X 10"^ 


0.81(30) 
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0.27(14) 


8 


1.1 


S3:F3 


0.091 < an < 0.15 


0.06(12) X 10"^ 


0.67(24) 


0.50(33) X 10"^ 


0.33(12) 


9 


1.0 


S3:F4 


0.091 < an < 0.16 


0.10(10) X 10"^ 


0.59(17) 


0.57(27) X 10"^ 


0.364(92) 


10 


0.95 


S3:F5 


0.091 <afl< 0.18 


0.54(73) X 10"" 


0.667(98) 


0.48(22) X 10~^ 


0.325(50) 


11 


0.87 



Table 6. Eq. (3.8) is fitted to the mode number per unit volume for sets SI and S3. Tlie fit range for set 
SI (amo = —1.15) is systematically explored. In fits S1:F1 to S1;F7, the fit-range higher end is kept fixed 
at 0.18 while the lower end is varied. In fits S1:F8 to S1:F26, the lower end is kept fixed at 0.091 while the 
higher end varied. S1:F4 (in bold) has been chosen as the final result. Fits S3:F1 to S3:F5 correspond to 
the lighter mass amo = —1.18. 7* obtained from this set is compatible with the heavier-mass result. 
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